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ABSTRACT. Much attention has been given to the problem of estimating 
cosmological parameters from the C; measured by future experiments. Many 
of the approaches which are being used either invoke poorly controlled ap- 
proximations or are computationally expensive. We derive exact results as 
well as fast and highly accurate approximations for mapping a theoretical 
model onto the observed power spectrum coefficients and computing their 
statistical properties. These results obtain from an analytic framework which 
applies for any azimuthally symmetric sky coverage regardless of the fraction 
of the sky observed by the experiment. 



1 Introduction 

There has been much recent at- 
tention to the problem of estimating 
the power spectrum of anisotropics 
in the cosmic microwave background, 
Ci, from observations (Tegmark 
1997; Bond, Jaffe & Knox 1998a,b). 
This is particularly topical because 
present and future experiments, from 
the ground, from balloons and from 
space, promise to provide a wealth 
of cosmological information (Knox 
1995; Jungman et al. 1995;Bond, Ef- 
stathiou, & Tegmark 1998; Zaldar- 
riaga, Spergel, & Seljak 1997). Even 
if we leave the formidable task of 
going from a time-ordered data set 
to a pixelised map to one side and 
assume the associated tasks of con- 
trolling systematic effects and fore- 
grounds will be so well controllable 
that they can be forgotten about in 
the final maps, the C; estimation 
still poses several fundamental diffi- 
culties. These stem mainly from the 
following facts: 

1. Balloon-borne and ground-based 
experiments observe only a frac- 
tion of the sky, while the theoret- 
ical C; are global quantities which 
would need observation of the full 
sky. Only in the full sky limit are 
the C; statistically independent 
quantities. 

2. Even satellite experiments cannot 
observe the full CMB sky due to 
galactic obscuration. 

3. Any measurement will be noisy, 
with a possibly anisotropic noise 
pattern. 

4. The sheer size of future data 
sets makes the estimation prob- 



lem computationally very expen- 
sive typically scaling as the cube 
of the number of pixels (Bond, 
Jaffe, & Knox 1998a,b). 

5. Even for Gaussian theories the 
Ci are non-Gaussian because they 
are quadratic quantities. In fact 
they are the sums of the squares of 
correlated Gaussian random vari- 
ables. The underlying Gaussian 
quantities are the mode ampli- 
tudes a; m of the spherical har- 
monic expansion of the sky signal. 

It has been realised long ago that 
the exact, theoretically trivial way 
of solving the problem by maximis- 
ing the Gaussian likelihood written 
in terms of the measured a; m is com- 
putationally not feasible. 

Therefore various approaches have 
been used to deal with these difficul- 
ties: 

Compress the data set in a more or 
less lossy fashion and then per- 
form a full likelihood analysis in 
the resulting system with fewer 
degrees of freedom. 

Neglect items 1 and 2 above , i.e. 
treat the C; as statistically inde- 
pendent. This is computationally 
a great simplification and hence 
very useful but leads to biased re- 
sults. 

Neglect item 5 above , i.e. treat the 
C; as Gaussian variates whose 
first 2 moments are identical with 
the true distribution. Thanks to 
the Central Limit Theorem, this 
is a good approximation as long 
as one can firmly establish the l- 
regime where it holds to the de- 
sired level of accuracy. The higher 
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order moments must decay suf- 
ficiently that the mode and the 
mean of the true likelihoods are 
coincident to a good approxima- 
tion. 

It should be noted that Oh, 
Spergel, &: Hinshaw (1998) have 
demonstrated an algorithm for solv- 
ing the likelihood problem efficiently 
by optimising their calculations ex- 
plicitly for the MAP satellite. This 
achieves a computational cost which 
scales with the square of the number 
of pixels. They also treat for the first 
time the masking of a small percent- 
age of the CMB sky by point sources. 
However, we believe that it is fair to 
say that their methods rely on the 
large sky coverage of a space mission 
for good performance and are still far 
from computationally trivial. 

The work presented in this talk 
presents a fresh look at the prob- 
lem. The goal is to step on to the 
path towards providing a mathemat- 
ical and computational tool which 
satisfies both the criterion of accu- 
racy to the exacting standards of to- 
morrows CMB experiments and is 
computationally feasible for a large 
and interesting class of cosmologi- 
cal theories, experimental setups and 
observational strategies. 



with 21 + 1 degrees of freedom and 
mean C theory . 

Owing to Galactic foregrounds, 
limited surveying time or other con- 
straints inherent in the experimen- 
tal setup, the temperature map that 
comes out of an actual measure- 
ment will be incomplete. In addi- 
tion, a given scanning strategy will 
produce a noise template. We model 
the noise as a Gaussian field Tjy 
with zero mean which is independent 
from pixel to pixel and modulated 
by a spatially varying rms amplitude 
Wfj (7) . Therefore the observed tem- 
perature anisotropy map is in fact 

f (7) = W(j) [T( 7 ) + W N (j)T N (j)] 

where W is unity in the observed re- 
gion and zero elsewhere. 

Expanding f as in Eq. (|l|) pro- 
duces a set of correlated Gaussian 
variates a; m for the signal and ajv i m 
for the noise. These combine into 
power spectrum coefficients 

Ci = ^ \aim + a N lm\ 2 (2) 

m 

whose statistical properties differ 
from the ones of the C;. Hence we 
refer to the C; as pseudo-Ci . 



2 Notation 

The full sky of CMB tempera- 
ture fluctuations can be expanded in 
spherical harmonics, Yj m , as£] 

T (l)= J2I2 a l™Ylm(f) (1) 

(=0 m 

where 7 denotes a unit vector point- 
ing at polar angle and azimuth 
<j>. A Gaussian cosmological theory 
states that the a; m are Gaussian dis- 
tributed with zero mean and spec- 

n . ^theory / 2\ 

lfied variance G ; = (|ci; m | ). 

Hence, for noiseless, full sky mea- 
surements, each measured C'i inde- 
pendently follows a x 2 ~ distribution 

1 We assume that there is insignif- 
icant signal power in modes with 
I > lmax and use the convention that 
sums over m run from —l m ax to l m ax 
and all quantities with index Im van- 
ish for m > I. 



3 Results: C; Statistics 

The approach we take is to work 
out the complete statistical predic- 
tions of Gaussian theories of struc- 
ture formation for the results ob- 
tained by a given experiment. The 
results we list in the following hold 
exactly for the case of azimuthally 
symmetric sky coverage (such as 
rings, annuli, polar caps, or a com- 
bination of the above, and "full" sky 
with a galactic cut) and a white 
noise distribution which is allowed 
to be anisotropic but must have the 
same axis of symmetry as the ob- 
served sky region. However, we show 
by comparison with extensive Monte 
Carlo simulations that our results 
are miraculously accurate even for 
strongly misaligned, symmetry vio- 
lating noise patterns, far into the 
noise dominated regime (so in this 
case read "approximate" in place of 
"exact" ) . 

We emphasize that while our 
Monte-Carlo simulations are com- 
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Figure 1. We show the theoretical standard cold dark matter C\ spectrum 
(solid grey line) superposed on 20 of our Monte Carlo simulations (dots). 
The white noise spectrum (grey short-dashed line) grows as i(/ + l)Cjv in our 
units. The long-dashed black line shows the sum of signal and noise while 
the white solid line displays the (pi) from 3328 Monte Carlo simulations. 
Apart from the large error due to non-Gaussianity at small I < 30, a small 
discrepancy can be seen for all I even for this COBE-style, almost full sky 
observation (ci also Figure^). 



puted for a ±20° galactic cut, our 
approach is not dependent on large 
sky coverage. We chose to compute 
this case as a demonstration of the 
accuracy of our formalism, its re- 
markable robustness with respect to 
a tilted noise template and its com- 
putational feasibility for maps with 
several million pixels. For such maps 
this is, to our knowledge, the only 
rigourous method which allows the 
discussion of such parameter biases 
and the computation of Ci statis- 
tics to the levels of accuracy which 
are quoted as the baseline for e.g. 
the Planck mission (Bersanelli et al. 
1996). 

Here are the results one by one 
(Wandelt, Hivon, & Gorski 1998a,b): 

1. An exact generalisation of the C; 
sample variance formula (Knox 
1995) to partial sky coverage and 
anisotropic noise. 

2. An exact analytical closed form 
solution for the non-Gaussian 
marginal probability distributions 
of the C;. Examples of these dis- 
tributions are shown in Figure ^| 

3. Exact expressions for cumulants 
and joint cumulants (and hence 
moments) of all orders for the C; . 
In Figure |s| we compare the true 



means, standard deviation, skew- 
ness and kurtosis to the ones for 
the x 2 distribution for all I. These 
can form the basis of a debiasing 
scheme, as suggested in Wandelt 
Hivon & Gorski (1998a). 

A new, standardised, x 2 dis- 
tributed statistic yi which mea- 
sures the goodness of fit between 
a set of observed C'i and a given 
theorjj^]. 

The discovery and removal of a 
correlation bias in cosmological 
parameter estimation which re- 
sults from neglecting the corre- 
lations induced by partial sky 
coverage and anisotropic noise. 
This is an effect which is dis- 
tinct from and in addition to 
the bias which comes from ne- 
glecting the non-Gaussianity of 
the Ci . This non-Gaussianity bias 
has been termed "cosmic" bias 
elsewhere (Bond, Jaffe, & Knox 
1998b). The correlation bias oc- 
curs at all I even those where the 
Central Limit Theorem ensures 
Gaussianity in the case of almost 



2 We thank G. Efstathiou for sug- 
gesting computation of this quantity 
to us. 
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Figure 3. Moment discrepancies between the WHG and x 2 distributions. 
We show the discrepancy in percent between the means, standard deviations, 
skewness and kurtosis for the Ci and \ 2 distributions. The stars are the (pi) 
computed from 3328 Monte Carlo simulations, showing excellent agreement 
only limited by Monte Carlo noise. Notice the oscillatory features in the 
discrepancy which indicate that it is dependent on the underlying theory, 
e.g. through the position and amplitude of the peaks. The difference in the 
means for I > 30 is dominated by the correlation bias and does not decay for 
large I. 



full sky coverage. ^|. To show the 
effect on the estimation of cosmo- 
logical parameters we form an ap- 



3 Of course our calculations also ac- 
count for the non-Gaussianity of the 
Ci distributions and are therefore 
free from "cosmic" bias. 



proximate likelihood by multiply- 
ing the marginal Ci -distributions 
and use this to estimate Qi,, the 
baryon fraction for the standard 
cold dark matter scenario. The re- 
sults are shown in Figure ^. While 
this is quite a small effect for 
the almost full sky case consid- 
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Figure 4. Correlation bias in cos- 
mological parameter estimation. We 
estimate the baryon fraction in the 
universe from 17 realisations of 
Ci, comparing the Gaussian and 
X 2 approximation to our likelihood 
(WHG). The true value is marked 
with the long solid horizontal line. 
The mean estimates are marked with 
short horizontal lines. Because we 
enforce the correct marginal distri- 
butions our result is guaranteed to be 
unbiased while the Gaussian and \ 2 
approximations neglect I to I corre- 
lations and hence produce a system- 
atic bias, here detected at more than 
3 sigma (dashed horizontal lines). 



ered in Wandelt, Hivon & Gorski 
(1998a) we conjecture this to be- 
come more important for medium 
to small sky coverage. 

The computation of any of these 
quantities, except for the joint mo- 

3 

ments, scales strictly as N^ ix , where 
N p i x is the number of pixels, with 
a small prefactor. This means all 
computation take on the order of 1 
minute for an i-range of to 1024. 
In the case of the joint moment 
yACjAC;/^ the computation scales 

as Np ix and takes on the order of a 
few minutes. 

Work is in progress to extend these 
methods to a viable tool for the 
joint estimation of the Big Bang 
Ci and cosmological parameters in 
Gaussian theories from noisy, cut- 
sky data (Wandelt, Hivon & Gorski 
1998b). 



This work was funded by the 
Dansk Grundforskningsfond through 
its funding for TAG 



References 

Bersanelli M. et al., 

1996, "COBRAS/SAMBA: Report 
on the Phase A Study", see also 



http: //astro, estex. esa.nl/Planck/ 



Bond J. R., Efstathiou G., Tegmark 
M., 1998, MNRAS, 291, L33 

Bond J. R., Jaffe A., Knox L., 1998a, 
Phys. Rev. D, 57, 2117 

Bond J. R., Jaffe A., Knox L., 
1998b, Radical compression of cos- 
mic microwave background data, 
[a.stro-ph/9808264 

Jungman G., Kamionkowski M., 
Kosowsky A., Spergel D. N., 1995, 
Phys. Rev. D, 54, 1332 

Knox L., 1995, Phys. Rev. D, 52, 
4307 

Oh S. P., Spergel D. N., Hinshaw 
G., 1998, An efficient technique to 
determine the power spectrum from 
cosmic microwave backk ground sky 
maps, astro-ph/980533£ , Ap. J. in 
press 

Tegmark M., 1997, Phys. Rev. D, 55, 
5895 

Wandelt B. D., Hivon E., Gorski 
K., 1998a, Cosmic microwave back- 
ground anisotropy power spectrum 
stati stics for high pre cision cosmol- 
ogy, |astro-ph/9808292| 

Wandelt B. D., Hivon E., Gorski 
K. M., 1998b, in preparation 

Zaldarriaga M., Spergel D., Scljak 
U., 1997, Ap. J., 488, 1 



Acknowledgments 



We would like to thank A. J. 
Banday for stimulating discussions. 



